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1. Biomass Properties and Prediction Tools for 
Vegetation Wildland-Urban Interface (WUI) Fires 


1.1 Overview 


Wildland fires are natural calamities that bring enormous environmental and economic damage worldwide 
and some of them cause human death. Achieving effective fire fighting is associated with the possibilities 
of predicting the characteristics of fire behaviour. Special attention needs to be drawn to meteorological 
conditions and their effect on fire spread in a bed of vegetation. On the other hand, studying the effect of 
forest fires on the environment is of equal significance for building the strategy and specific actions to be 
undertaken in the fire fighting process. Especially, Mediterranean climate countries are subjected to higher 
risks of fires and all the damage they involve. 


This material presents a review on related literature involving main characteristics and fire behaviour 
prediction models for surface fire and especially for pine litter species. The first part of this section 
presents experimental data for various characteristics of pine needle species, and also results from 
laboratory observations and studies for fire behaviour in a vegetation layer comprising above species. The 
second part gives consideration to various types of models used to predict the behaviour of surface fires. 
Also presented are basic approaches employed in describing the processes involved in the heating up, 
ignition and burning of the vegetation, and also the effect of these processes on the parameters of the flow 
in the fire zone. Data for model verification are also presented. 


1.2 Research in the Field of Occurrence and Spreading of Forest Fires 


Three components should be simultaneously available in order for a forest fire to occur: fuel (forest 
vegetation), oxygen and a heat source. 


Fuels in the forest can be divided into four heights (Missbach et.al., 1982): 


e Vegetation of a height of above 2 meters, 

e Shrubs and low-height trees up to 2 meters high that are most commonly encountered in the 
Mediterranean region. Major representatives of these species are various shrubs and grasslands; 

e A layer of dry grass and tree leaves litter; 

e Soil cover bed. 


Fires occurring in the last three groups are considered surface fires. 
Experiments carried out to examine the vegetation’s susceptibility to ignition and the behaviour of the 


different types of vegetation as a fuel bed resulted in final output and obtained values for the so called fuel 
particle characteristics. The latter are among the basic data used to design the fuel models. 
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On the other hand, the studies and the assessment of the intensity and the risk of occurrence of forest fires 
are based on the determination of fire behaviour characteristics. The most common task involved in 
obtaining the rate of spread of the fire under various surrounding conditions and for different types of 
vegetation species is to establish the current location of the fire outlines along a given direction and also 
the average rate of variation of these outlines — the rate of fire spread (ROS). 


In order to combat successfully the ignition of forest vegetation and to fight the fires that occur, one of the 
conditions is the availability of adequate tools to predict the behaviour of forest fires. Various models are 
being developed and used to act as tools to determine individual fire behaviour characteristics or 
parameters in the fire zone. These models are obtained using various approaches and these vary from 
statistical processing of experimental data to the presentation of a detailed physical and mathematical 
picture of processes associated with the occurrence and spreading of the fire. 


1.2.1 Strengths and Weaknesses of Literature 


Generally, a literature survey of the problematic of a study allows to acquire information which can then 
be used as a basis for analysis of problems and to outline the areas where the main efforts and resources 
should be concentrated. 


A large amount of information is available in the readily accessible literature sources for fuel particle 
characteristics (mainly physical ones) for different vegetation species of low and medium height, 
characteristic of the Mediterranean European regions. This also provides a possibility to determine these 
characteristics by means of studies using samples in laboratory conditions and employing simple methods 
and settings. 


Another strength of the literature sources is the collected information and the studies carried out with the 
aim of developing models to be used for predicting the outlines of the fire and the fire rate of spread based 
on statistical processing of data on these characteristics. The most intensive studies in this area and, 
respectively, the most comprehensive data is to be found for fires in vegetation beds comprising shrubs 
from the region of Australia, as well as some studies for North America — the US and Canada. 


The possibility for using a data base which compiles numerous reports from studies that have been carried 
out is also one strength of the literature survey. The EUFIRELAB data base provides comprehensive 
information on activities involved in the prevention of forest fires in the Mediterranean countries. This 
data is structured observing the relevant characteristic attribute feature and provides possibility to trace the 
developments in the study and the achievements in the relevant areas. Another advantage is the possibility 
to establish direct contact with the authors. Another powerful data source is the International Journal of 
Wildland Fire and this is the formal scientific edition of the International Association of Wildland Fire. 
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Weaknesses involved in the available literature sources for experimental studies lie in several basic 
aspects. Data on fuel particle characteristics for each individual type of vegetation species vary (in terms 
of number) within wide ranges. This is true for both data provided for the same region, and also for 
comparison studies for several different regions. These differences are due to the variations in the 
characteristics of vegetation particles resulting from the different climatic conditions, as well as to the 
different methods and methodologies adopted for determining those characteristics. Data on the methods 
employed is just as scarce and hinders the assessment on the applicability of these literature data. 


Determining fire behaviour characteristics under laboratory conditions is suitable for small scale fires 
which comprise the group of fires occurring in low height (depth) vegetation beds. Therefore, laboratory 
experiments are only limited to studying mainly fires occurring in pine litter other then forest litters. A 
limitation in these experiments is also the space (time) for modelling the phenomena and hence, for 
obtaining data for large scale fires and well established fire behaviour characteristics. The experiments for 
modelling surface fires into “live” vegetation layers (shrubs, grasslands) need to be performed in actual 
conditions and this makes it difficult to organize them and to provide the necessary measurement 
equipment, as well as making them more expensive. Therefore, literature sources providing evidence for 
such type of experiments are very rare. 


Data on thermal decomposition and combustion of different vegetation species is also very scarce in the 
literature sources considered herein under. Experimental studies in this area are hindered by the need for 
using complicated and expensive laboratory equipment and installations, as well as by the complex nature 
of the processes involved in the modelling of close to actual conditions. The theoretical modelling of such 
processes is rather complicated and they need to be simplified in order to obtain results that are applicable 
for the relevant field. 


Studies in the field of fire behaviour for vegetation occurring in the region of Greece are very limited and 
much less for other Mediterranean European countries, such as Italy, France, Spain and Portugal. 


According to Allgower et al. (2002) some of the gaps in modelling the fire behaviour are linked to the fact 
that ignition and maintenance of the fire are a function of various conditions determined in turn by an 
overall fuel complex. However, in experimental studies the attention is concentrated basically on the 
formation of fuel models to be used to describe the relevant type of vegetation species. However, the latter 
are only part of the overall fuel complex. Another drawback of the modelling of forest fires is the fact that 
propagation modelling is not yet really ready for 3D-fire spread modelling over vast areas. This is 
preconditioned both by insufficient theoretical studies in this field, and also by the high computational 
cost pf the numerical experiments and this becomes an obstacle for shifting from small-scale combustion 
modelling to larger-scale 3D-fire propagation modelling. Another gap encountered in fuel description and 
modelling is how surface fires spread into active crown fires. 


1.2.2 European Programs 


A major role in the development of research programmes in the area of forest fires for the region of the 
Mediterranean countries is played by the organization EUFIRELAB (www.eufirelab.com). The 
EUFIRELAB is a Euro-Mediterranean wildland fire laboratory, a wall-less laboratory for wildland fire 
sciences and technologies in the Euro-Mediterranean region. 


Download free ebooks at bookboon.com 


Biomass Properties and Prediction Tools for 
Biomass Properties and Fire Prediction Tools Vegetation Wildland-Urban Interface (WUI) Fires 


EUFIRELAB is structured in units aimed at reinforcing the co-operation among Euro-Mediterranean 
teams, activating large exchanges of knowledge and know-how, developing common concepts, 
approaches, and “languages” and fostering the common use of facilities for research and/or technological 
development. The major fields of research on the issues of forest fires (as indicated with codes), 
distributed among individual research teams are: 


1. WP02 Wildland fuel description and modelling unit 
e WP02T1: State of the art and survey 
e WP02T2: To elaborate common methodologies 


2. WP03 Wildland fire behaviour modelling unit 
e WPO3TI: State of the art and survey 
e WP03T2: To compare the different types of models 
e WPO03T3: To define common procedures for 
e WPO03T4: Towards a European scale of magnitude of wildland fires for characterising the 
intensity of wildland fires and prescribed burnings 


3. WP04 Wildland fire, ecosystems functioning, and bio-diversity unit 
e WPO04TI: State of the art and survey 
e WP04T2: Fire impacts on the different components of the ecosystems 
e WP04T3: Methodologies and tools for analysing and monitoring vegetation dynamics and 
restoring burned areas 
e WPO04T4: Prescribed burning, a tool for managing bio-diversity and ecosystems functioning 


4. WP05 Socio-economy unit 
e WPOSTI: State of the art and survey 
e WPO05T2: Towards common methodologies for investigating on wildland fires causes and factors 
e WPO0S5T3: Towards common methodologies for studying the costs-to-benefits ratio of wildland 
fires prevention 


5. WP06 Decision support tools unit 
e WPO6T1: State of the art and survey 
e WP06T2: To elaborate common specifications for further decision support systems 


6. WP07 Widland fire metrology unit 
e WPO7TI: State of the art and survey 
e WP07T2: Towards common methodologies for collecting data during laboratory experimental 
fires 
e WPO07T3: Towards common methodologies for collecting data during outside fires 
e WPO07T4: Centre for technological development 


7. WP08 Widland fire risks and hazards unit 
e WPO8TI: State of the art and survey 
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e WP08T2: Common methods for mapping wildland fire risks 
e WPO08T3: Towards a Euro-Mediterranean Wildland Fire Danger Rating System 


8. WP09 Wildland fire suppression unit 
e WPO9TI: State of the art and survey 
e WPO09T2: Towards common methodologies for wildland fire suppression planning 
e WPO9T3: To improve the safety and efficiency of the fire fighters 


10. WP10 Widland — urban interfaces management unit 
e WPIOTI: State of the art and survey 
e WP10T2: Towards common methodologies for managing wildland-urban interfaces 


Morvan (2004) present the basic characteristics of the initial versions of ten codes of behaviour models, 
developed in cooperation with EUFIRELAB - FIRESTAR-2.0, FireRegime-1.0, SPREAD Section2, 
SPREAD Section2 WP, FireStation, FIRE LINE ROTATION MODEL (FRM), INCENDIU 1.0, 
SPREAD 1.0, AIRFIRE and DISPERFIRE. 


One of the projects financed by the European Union is EFAISTOS. It constitutes the first part of the 
European behaviour model for forest fire, taking into account the lower layers of the forest fuels. Present 
In Valette et al. (1999) is presented a generic model of the experimental process, based on expert 
knowledge and experiments carried out within the framework of the European R&D project "Efaistos" 
(Environment and Climate Program, DG XII). This model can be used in some experimental fires and 
some real fires (prescribed burning and wildfires). It has been implemented within a Fire Behaviour 
Modelling Support System, aiming to constitute a database at the European scale. In order to assist the 
researchers in their activities related to the experiments, the researchers propose to design a prototype of 
Problem Solving Environment (PSE), called MODELLIS (MODELLing Information System), dedicated 
to the fire behaviour modelling support founded on terrain experiments. 


PROMETHEUS is a Management Information System dedicated to forest fire protection. The project is 
realised in the frame of ENVIRONMENT and CLIMATE R&D Program of the European Union and was 
implemented by a Consortium of seven European research organisations co-ordinated by 
ALGOSYSTEMS s.A 


1.3 Fuel Particle Characteristics 


1.3.1 Description of Fuel Properties 


Vegetation is assimilated to fuels for the purpose of predicting the behaviour of a fire, either relative (fire 
danger rating) or actual (rate of spread, intensity, flame size, etc.), as well as the impact of fire. 
Consequently, there exists an overwhelming variety of methods designed to study and describe fuel 
characteristics. 
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Fuel particles are the smallest elements considered in order to study the fuel structure. They are organs or 
pieces of the aerial parts (dead or live) of the vegetation: branches, leaves, barks, cones, needles, etc. Fuel 
particles are partitioned by size classes and their condition (live or dead), therefore establishing the limits for 
the description of their properties. Particle size categories are not standard across the world, which hinders 
comparability of fuel data and fire behaviour models. For fine fuels the limit value for diameter is 6mm. 


Fuel particle characteristics contribute to the prediction of wildland fire intensity and severity, with all its 
consequences on suppression difficulty and human safety. Characterization of fuel particles is therefore 
required to interpret the results of flammability experiments in the laboratory and as an input to 
semiempirical and physical fire behaviour models. 


Main types of characteristics of fuel particles are: 
1. Physical characteristics: 
- measured physical characteristics: length, width, thickness, diameter, mass, volume; 
- calculated physical characteristics: mass to volume and surface to volume ratios. 
2. Chemical characteristics: moisture and chemical composition, ash content. 


3. Thermal characteristics: thermal degradation and high calorific value. 
4. Other: flammability. 
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The physical, chemical and thermal properties of fuel particles are assessed at the level of the individual 
particle or element (leaf, spine, stalk, twig, branch, stem, etc.), or of compounded particles belonging to 
the same biological entity, e.g. the assemblage of leaves and small twigs of a given shrub species. 


Extensive search in biomass has been also made in order to assess the relative flammability of the species. 
Some authors use the term ‘inflammability’ as the ability of the fuel to ignite after having been submitted 
to calorific energy. This term coincides with the term ‘ignitability’ in the American literature. According 
to fuel flammability the species are ranked by using two properties — heat content and temperature of 
ignition (Dimitrakopoulos, 2001). 


The following fuel properties affect the plant flammability: 


Heat content (calorific value) of a given fuel is a comprehensible measure of the potential thermal 
energy that can be released during the burning of the fuel; 


e Total ash contents reduces the amount of combustible fuel mass, since only the organic part of the 
fuel supports combustion. Mineral (silica-free) ash is important during the early stages of 
pyrolysis by catalyzing the char formation reactions, thus reducing the flammable gases evolved. 


e Surface area-to-volume ration is a meaningful measure of the fuel particle size since it determines 
to a great extent the heat and moisture exchange rates. 


e Fuel particle density affects the thermal conductivity of the fuel and, therefore, its ignition time. 


Various methods and procedures are available for sampling and for determining fuel characteristics. 
Presented underneath is an example based on data derived from L. Nunez-Regueira (1996). 


For collection and preparation of the samples, 1 ha of woodland is chosen. The plots are divided into 1 m? 
size sites, five of which are randomly chosen. From every site, bulk samples consisting of bark, branches 
having a diameter not greater than 8 cm, fruits, leaves and in general all of the living parts of trees are 
collected. The bulk sample is reduces by coning and quartering procedure to a representative sample of 
about | kg. Part of the sample is used in the flammability experiments. The environment of the fuel 
samples is recorded and example is presented in Table 1.1. 
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Item Data 
Altitude 100-320 m 
Annual rainfall index 1865 mm 
Summer rainfall index 217mm 
Mean annual temperature 12.1 °C 
Mean daily maximum temperature of the warmest month 315°C 
(June) 
Hydric deficiency 26 
Mediterraneanity index 1.27 
Representative species of the zone Pinus Pinaster Aiton, 


Eucalyptus globulus Labill , 
Sarothamnus scoparius (L.), 
Ulex europeas L. , 

Rubus fructosus L., 
Pteridium aquilinum L., 
Castanea sativa Miller, 
Quercus robur L., 

Acer pseudoplatanus L. 


Table 1.1 Environment of the fuel samples. 


1.3.2 physical, chemical and thermal characteristics 
1.3.2.1 Methods for determining the vegetation particle characteristics 


In their publication Allgower et all. (2002) present the basic methods for experimental determination of 
various characteristics of the vegetation bed in terms of using it as a fuel bed. The following presentation 
summarises these methods and is based on this report of EUFIRELAB. 


The density of fuel particle is usually assessed by weighing ovendry particles and by measuring the volume 
of fresh, air-dry or oven-dry particles, although some researchers use fresh samples. Volumes are measured 
by immersion in mercury or in water, the most common choice. There is also another approach where the 
density is determined based on measurements of particle size and calculations of particle volume. 


Download free ebooks at bookboon.com 


12 


Please click the advert 


Biomass Properties and Prediction Tools for 
Biomass Properties and Fire Prediction Tools Vegetation Wildland-Urban Interface (WUI) Fires 


Determination of the surface area to volume ratio of fuel particles is possible by a variety of methods. Like 
with particle density, measurements can proceed in fresh, air-dry or oven-dry biomass. The most accurate 
solution for three-dimensional and long narrow particles is the measurement of the perimeter-area ratio of 
a cross section by photomicrography or image analysis. The most straightforward and used approach is to 
establish mathematical relationships between the physical dimensions of fuels and their surface and 
volume which rely on the description of shape by simple geometry Johnson (1984) propose a method for 
conifer needles based on simple geometry where the surface area is a function of needle length, number of 
needles per fascicle and volume displaced by water immersion 


Various methods could be used to determine the moisture content of the fuel bed. The most commonly 
used one among them is based on measurements of vegetation mass before and after the drying process, 
and also on calculations of moisture content as the difference between the two values and then expressed 
relative to the mass of either the moist vegetation or the dried vegetation. The use of portable ovens, either 
based on microwave radiation or on the conventional method, can overcome the major problem of oven- 
drying, i.e. the time delay to obtain a result and its impracticability in the field. Chemical methods for 
determining the moisture content of the vegetation bed were also established and these involve, for 
example, the addition of calcium carbide to minced fuel samples in a pressure cylinder. The moisture 
content of the fuel bed can also be determined by measuring the electrical characteristics of the vegetation 
particle (capacity or electrical resistance) and using the relation between these quantities and the amount 
of moisture content. 
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Determination of ash content can be performed with the samples placed in a muffle furnace heated to 
temperatures higher than 300°C and up to 650°C, cooled down and then weighed. An alternative method 
for determining this chemical characteristic is to use a differential thermal and thermogravimetric analysis 
which involves the heating up of the fuel particle with a constant rate, and then measuring the temperature 
and the variations in the mass of the sample. These methods can be used for determining also the chemical 
composition (the amount of volatiles, etc.), and also to determine the rate of variation of the fuel mass. 
Heat contents of fuel particles, whether high or low (i.e. corrected for the latent heat of water vaporisation) 
and with or without ash, are determined by standard adiabatic bomb calorimetry. The only alternative 
method has been developed by GILLON et al. (1997) and uses near-infrared reflectance spectroscopy. 


1.3.2.2 Experimental data for Mediterranean pine species 


Presented further are data for the characteristics of pine needles for vegetation species typical of the 
Mediterranean region. Most of the data is characteristic of a fuel bed made up of mainly dead vegetation 
particles. Table 1.2 provides data from experiments carried out using EUFIRELAB programmes for the 
countries along the Mediterranean for three major representatives of pine species — Pinus halepensis, 
Pinus pinaster and Pinus brutia. 
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Table 1.2 Characteristic properties of pine species from different sources (in parentheses) 
* surface area to volume ratio, [1/m] 


Pinus brutia and Pinus halepensis have high heat content, surface-to-volume ratio and very low ash 
content and particle density. Since plant flammability is favorably affected by SA/V and heat content the 
pointed species (Pinus brutia and Pinus halepensis) are considered very flammable (Dimitrakopolous, 
2001). Pine leaves have smaller SA/V ratio then needles and greater mass-to-volume ratio (density). It is 
obvious from the results presented that variations in heat and ash content values for different sources are 
not substantial, while substantial differences are observed for surface-to-volume ratio and density. The 
latter could be due to different factors: the particular region where the study was performed, the season of 
the year, the moisture content of vegetation, etc. Small differences in heat and ash content could also be 
interpreted as a sign of very similar chemical composition of these vegetation species. 


The tables underneath present data from studies on variations of fuel particle characteristics for Pinus 
pinaster (Nunez-Requerira, 1996). The studies have been conducted in the region of La Coruha (Galicia, 
Spain), for two different areas - Sada (coastal area) and Santiago (hillside and plateau area). These data 
were evaluated as a help for fighting forest fires, which have been very frequent in this region. Table 1.3 
illustrates the variations in fuel moisture content (expressed as a percentage of the initial mass of 
vegetation particles), density and ash content. As it is rather normal to expect, the lowest moisture content 
is observed in sprint time, but high values (above 50%) should also be noted, most probably due to the fact 
that live foliage is studied and the area is very close to the Ocean. 
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Season Moisture, Density, Ash, 
% kg/m? [%] 
Spring 56.30 650 1.04 
Summer 53.50 640 0.84 
Autumn 61.50 640 0.91 
Winter 60.46 630 0.71 


Table 1.3 Pine litter fuel properties at different seasons — Pinus pinaster needles. 


When presenting the results on heat content variations throughout the year two calorific values are 
considered. The higher heating value (HHV) is defined as the quantity of heat generated by complete 
combustion, in a bomb calorimeter, of a unit mass of a sample in an oxygen atmosphere assuming that 
both the water contained in the sample and that generated from the combined hydrogen remains in liquid 
form. The lower heating value (LHV) can be calculated if it is assumed that the water in the products 
remains in the form of steam. Values for HHV and LHV are indicated in Table 1.4. 


Season HHV, LHV, 
[kJ/kg] [kJ/kg] 
Spring 19480.91 6524.90 
Summer 20658.80 7646.12 
Autumn 20463.05 5868.74 
Winter 20398.48 6028.93 


Table 1.4. Heating value of fuel as function of the season. 


As expected these calorific values are highest in summer and lowest in spring. As was mentioned, this fact 
is related to the blooming period of most of the species, that Mediterranean Greece coincide with a season 
of frequent rain fall, so increasing their moisture content and diminishing their HHV. The experiments 
show that these values differ along the place of the probe, which is related with the variety of other species, 
i.e. trees and bushes with high calorific value in the appropriate zone. As can be seen Pinus pinaster has 
extremely high calorific power and also high flammability during the whole year. The great amount of 
essential oils and resins, combined with a HHV and the high inflammability of this species, places it in the 
group of high-risk trees because of its readiness to start and spread forest fires. 


Wood can be chemically analysed by breaking it down into structural components (called proximate analysis) 
or into chemical elements (ultimate analysis). The main structural components are Cellulose, Hemicellulose, 
Lignin and ash. Volatiles and solid carbon (char) are generated in the process of heating-up the fuel. Carbon 
and oxygen are the main chemical elements constituting almost 90% of the vegetation contents. 
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Nunez-Requerira (1996) also presents data on the chemical composition of Pinus pinaster. The chemical 
analysis (ultimate analysis) of vegetation species is presented in Table 1.5. 


Chemical analysis Volatile metals 
Season [% of total composition] [ppm] 
N Cc H O Ss Cl Cu Cd Zn Pb Mn 


Summer | 2.25 | 50.13 | 6.54 | 40.46 |0.10 |0.52 | 52.34 1.94 |19.40 |6.78 | 6390.00 


Autumn 1.90 | 44.64 |5.71 | 47.46 |0.22 | 0.067 | 6.13 2.30 | 17.48 1.86 | 83.12 


Winter 1.65 | 48.10 | 7.31 | 42.83 |0.06 /0.05 | 4.32 2.88 | 11.53 0.90 | 92.24 


Spring 3.02 | 46.96 |6.39 | 43.41 |0.22 |- 0.29 1.60 | 25.60 2.10 | 1206.30 


Table 1.5 Ultimate analysis and volatile metal analysis of Pinus Pinaster. 


As can be observed from the table, the composition of the fuel does vary in small ranges during the year, 
but generally the chemical composition could be considered constant (Mathieu, 2002; Viegas, 2001). 
Furthermore the fuel shows its highest flammability in the summertime and lowest in spring and winter as 
expected. These last values relate to the high degree of humidity of the species in these seasons. The 
observed high contents of Mn, compared to other elements can be explained by the need for this ion in the 
transportation from the water to photosynthesis areas. 
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For comparison proximate and ultimate analysis of some other biomass fuels are presented in Table 1.6. 


a alana Proximate analysis 
[% of total composition] ¥ 
Fuel Moistur | Ash, ; 
Fixed |HHV | LHV 
Cc H |O N S Cl e, wt% | wt% dry , 
wet fuel | fuel ° De S| 
clean 49.10 | °° |443 | o.48 | 0.01 }0.10 | 50 1.3 13.2 |192 |15.4 
wood 0 0 
Bark fi 
ark from | 54 49 [89 | 42-4 | o 44 10.03 10.03 |55-65 |2.34 |2246/19.83 | 18.54 
spruce 4 0 
Straw from | 49 69 | ®1 | 43-5 | o61 |0.07 |0.18 | 55 471 17.59 118.94 | 17.65 
wheat 6 0 
2/427 
Iss TEeO ae ag, | 1.54 10.15 |0.07 |60 8.85 |17.65 |18.37 | 17.13 
canary 5 0 
Waste— 
2 | 45. 
demolition | 48.80 : 0.15 |0.03 }0.08 | 15-20 |0.9 : 154 | 13.9 
wood 
4 
Salix 50.30 : 43.1 | 0.40 | 0.03 | 0.004 | 50 1.18 |18.92 |19.75 | 18.42 


Table 1.6 Comparative table for the properties of different biomass fuels. 


Additional data for the chemical composition of the pine needles is presented by Ucar (2004). This data 
reveals the chemical complexity of the fuel, the needles, the branches as well as other components of the 
fuel, so the expected chemical reactions set will consists of hundreds of reactions, which are not 
applicable in an engineering approach to model a wild fire behaviour. Furthermore, the composition of the 
released light gases (volatile gases) as well as released tar varies with the environment conditions (Barnola, 
2000; Owen, 2002; Dormont, 1998; Piccardo, 2005) — temperature of the fuel particles, intensity of 
heating rates to/from the fuel, gas composition environment (mainly the partial pressure of the oxygen). 
Obviously the presented, strongly time dependent process can not be described in deep details in any of 
the available biomass combustion models, both because of the uncertainty of the initial and boundary 
conditions as well as the fuel properties, and the complexity and uncertainty in the chemical reactions and 
their rates. 


As can be seen the variety of biomass fuels presents close range of fuels properties variations, thus 
revealing opportunity to develop a common model for biomass combustion (Carvalho, 2002; Jones, 2000; 
Demirbas, 2004), based on the appropriate physical properties and chemical composition of the fuel 
(Owen, 2002; Dormont, 1998). This will be used to develop a model, based on data for different biomass 
materials in case of lack of specific information for items of the developed model. For example the HHV 
could be calculated, based on the ultimate analysis as follows: 
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HHV = a GHC +102H +6.3N +19.1S —9.850), [MJ/kg dry fuel] (1.1) 


Fletcher (2000) proposes a formula, which is used to calculate the HHV of the biomass fuel as follows: 
HHV = 2.32(146.58C + 568.78H —51.53(0 + N)—6.58Ash + 29.45), [kJ/kg], (1.2) 


where C, H, O, N and Ash — weight percentage on the dry basis. 


Another formula (Fagbemi, 2001) is proposed by Institute of Gas Technology (IGT), in which the amount 
of the elements, C, H, O, N and Ash, are expressed in mass fractions: 


HAV = 354.68C + 1376.29H + 71.26 —15.92 Ash —124.69(0 + N) [kJ/kg], (1.3) 


Actually the HHV for pine litter fuel is available and the presented formulas could be used for validation 
and verification purposes as well. 


1.3.3 Ignition Temperature and Flammability 


The second phase (after pyrolysis) in the burning process is exothermic and it is known as combustion. 
Combustion may be with or without flame. Flaming combustion is a gas oxidation phase, accompanied by 
emission of flames. It takes place when the temperature of volatiles, emerging through the surface of the 
forest fuels, reaches 450-500 °C. Liodakis (2002) presents data, according to which the minimum surface 
temperature (critical surface temperature) of wood under radiative heating mode for spontaneous ignition 
is 600 °C and for piloted ignition 300-410 °C. With convective heating the spontaneous ignition, occurred 
at 490 °C and with piloted ignition 450 °C (Dormont, 1998). From the literature is shown that both the 
minimum pyrolysis rate and the minimum surface temperature can be used as criteria for ignition. 
However, from a practical point of view, the surface temperature criterion is much more convenient. 


This same report also presents results from experiments for ignitability properties for Pinus halepensis and 
other vegetation species characteristic for the particular region of Greece. Form these we also observe that 
an ignition temperature of 480 °C was defined for Pinus halepensis. Moreover, the generation of volatiles 
begins at temperatures above 200°C and it reaches its highest intensity in the range of 320 — 370 °C. When 
a temperature of 500°C is reached, the generated amount of volatiles constitutes around 55% of the initial 
mass of the sample material. 


Flammability (the ability of species to ignite and sustain fire) ranking of forest species is an essential 
component of fuel hazard and fire risk assessment. Dimitrakopoulos (2001) presents some studies of 
above-mentioned vegetation species characteristic for the particular region in Greece. Pinus pinaster and 
Pinus brutia, which dominate the composition of pine litter, have been classified as “Extremely 
flammable”. Table 1.7 presents appropriate ratings according to the flammability index (L. Nunez- 
Regueira, 1996). 
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Index Flammability 
0 Very low flammability 
1 Low flammability 
2 Flammable 
3 Moderately flammable 
4 Very flammable 
5 Extremely flammable 


Table 1.7 Flammability index. 
Table 1.8 illustrates the seasonal variations in the flammability index of Pinus pinaster throughout the 
year (L. Nunez-Regueira, 1996). The expected highest fire risk is during the summer. Generally, from this 
data it can be seen that pine species are extremely flammable and fires occurring in such vegetation 
species are usually highly intensive and carry a high degree of risk. 


Flammability values of Pinus Pinaster 


Season Spring Summer Autumn Winter 


Index 4 5 4 3 


Table 1.8 Flammability of Pinus pinaster at different seasons. 


1.3.4 Development of Fuel Models — Basic Characteristics 


Fuel complexes result from the organization of fuel particles into a microstructure, which can result of 
single, or multiple beds or layers. 


Fuel group is a complex of fuel elements with average properties values representative of the typical fuel 
conditions (combustible materials) of a certain vegetation type (Dimitrakopoulos, 2002). It is also a set of 
quantitative fuel inputs to fire behaviour models. 


A major characteristic feature of the fuel group is fuel load - the dry weight of fuel per unit area, [kg/m’]. 
The structural arrangement of fuel particles within a fuel complex is essentially defined by porosity, 
particles orientation, and vertical distribution. Porosity is the ratio between the void volume in a fuel 
complex and its surface area. Fuel depth equals fuel height in grassy or shrubby fuel layers, but if the fuel 
complex includes a litter bed (usually not entirely available to burn actively in the fire front), a more 
objective definition of fuel depth is the vertical extent of the combustion zone. 


Particles orientation and vertical fuel distribution are not explicit factors in fire models. 


Download free ebooks at bookboon.com 


21 


Please click the advert 


Biomass Properties and Prediction Tools for 
Biomass Properties and Fire Prediction Tools Vegetation Wildland-Urban Interface (WUD) Fires 


According to the classification undertaken in Dimitrakopoulos (2002), 7 fuel models were defined for the 
region of Greece. Dominant species for fuel model 7 (forest litter layer) are Pinus Brutia and Pinus 
Halepensis. Other models are for medium to high shrubs and one model is provided for grassland. Fuel 
models 1 — 5 include the following vegetation species Pistacia lentiscus, Quecus coccifera, Arbutus unedo, 
Plomis fruticosa and Sarcopoterium spinosum. 


The methodology for the grouping of categorically equivalent species into fuel groups in Greece consisted 
of the following steps: 


(i.) All areas covered with Mediterranean vegetation in Greece were stratified on vegetation maps 
according to the dominant vegetation type: grasslands, phrygana (small, xeric shrubs up to 0.5 m 
height), maquis (evergreen-broadleaved, sclerophyllous shrubs, 0.5-3 m height), or closed-forest 
litter of Mediterranean pine species (Pinus brutia and Pinus halepensis), the latest of main interest 
in this work. 


(ii.) In every representative location, 12 fuel parameters are measured in 500 m’ sampling plots. The 
clip-and-weight method is used for the determination of all fuel loads by size category. The line- 
intercept method is used in order to estimate the area cover by each vegetation type. All fuel loads 
(fuel weight per unit surface area) are expressed on a dry-weight basis. 


(iii.) | The collected data are subject to statistical analysis. 
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Fuel categories for the region of Greece drawn up by Dimitrakopoulos (2002: 128) can be divided into 
three separate groups: 


e Shrubs of average height between | and 2 meters — the first three models; 

e Shrubs of average height between 0.4 — 0.55 meters — the Phrygana I and Phrygana II models 
(here we also include the grasslands model); 

e Litter representing models —Pine forest litter. 


At present we are interested to develop the conditions for the burning of Pine litter. The literature survey 
shows that the forest litter of Mediterranean pine species has the properties, presented in Table 1.9. The 
chemical composition of the pine needles will not be presented here, but can be seen in Ucar (2004). 


Forest litter of Mediterranean pine species 
Average _| Branch diameter, [cm] Litter Litter Dry pine ious 
; % cover tree- 
height, depth, weigh, SE needles, ppsuenes 
[cm] 0.0-0.5 | 0.6-2.5 | 2.6-7.0 [cm] [t/ha] y [t/ha] , 
[t/ha] 
1.5 70 % 20 % 10% 6 12.55 100 10.2 2.35 


Table 1.9 Properties of forest litter of Mediterranean pine species 


The distinct morphological differences of the two dominant phrygana species in Greece resulted in two 
separate fuel models with different geographic distribution. The low-elevation Mediterranean grasslands 
and the litter of closed pine forests demonstrated limited spatial heterogeneity and are represented by a 
fuel model each. The fuel model for pine forest litter considers total load of 12.55 t/ha. Due to the fact that 
the experiments for the fuel properties determination are held in the summer period, most of the grassland 
fuel load was allocated to the dry fine fuel category. The litter load of the Mediterranean closed forests of 
Pinus Halepensis and Pinus Brutia (Calabrian pine) comprised mainly dry pine needles (10.2 t/ha) with a 
small proportion allocated to the fallen tree-branches (2.35 t/ha) of the forest floor. Experiments show that 
severe burning conditions are observed in the forest litter despite the compactness of the litter layer (litter 
depth 6.0 cm). It should be emphasized that the local fuel conditions may vary from one area to another, 
thus the presented fuel model may substantially vary according to the fuel data. 


1.4 Theoretical Modeling of Fire in Vegetation Layer 


1.4.1 Classification of Fire Behaviour Models 


Various systems for classification of models used to predict surface fires’ behaviour are available. One of 
these classification systems groups the different models according to the predicted characteristics (Grishin, 
1997): 


e models used to predict the rate of spread (ROS) of the fire; 
e models used to predict the outlines of the forest fire; 
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e models used to predict the flow characteristics, the heat transfer and mass transformations along 
the fire front and inside the fire zone; 


e general mathematical models within which all characteristics could be predicted (rate of spread, 
outlines of forest fires, temperature range, component concentration and velocities distribution) 
along the fire front and inside the zone of the forest fire. 


Depending on the approach adopted in the development of the models used to describe the behaviour of 
forest fires, including that of surface fires, two basic groups of models were formed: stochastic models and 
deterministic models. Stochastic models consisting to predict the more probable fire behaviour from 
average conditions and accumulating acknowledges obtained from laboratory and outdoor experimental 
fires. Deterministic models in which the fire behaviour is deduced from the resolution of the physical 
conservation laws (mass, energy, momentum) governing the evolution of the system formed by the flame 
and its environment 


The main purpose of these models is to predict the local rate of spread of a fire front, when parameters 
characterising the condition of spread (vegetation, meteorology, terrain) are given. 


It is reminded that stochastic models are only based on the observation of field fires (experimental fires 
and wildfires) from which the fire rate of spread (ROS) is related to relevant parameters in a purely 
statistical way (fuel type, fuel loading, fuel moisture, wind). These empirical relations depend strongly 
from the very specific conditions from which the statistical study was performed. Without systematic 
parametric studies, it is very difficult to extract a general behaviour for the fire. 


The most commonly used classification of fire behaviour models groups the latter into three basic 
categories: empirical (or statistical), semi-empirical (semi-physical or laboratory models), and physical 
(theoretical or analytical) (Morvan et.al., 2003). The purpose of physical modelling of fire behaviour is to 
obtain a mathematical solution to the complex mechanism of the occurrence and spreading of forest fires 
(Andrews, 2001). Formulating an accurate theoretical model capable of providing a reliable description of 
fire behaviour that is of certain practical value requires multiple combinations of input conditions (to 
comply with actual ones) and powerful computing resources. 


At present, fire modelling for operational applications is basically addressed to the establishment of semi- 
empirical or empirical models. 


1.4.2 Empirical Models 


Empirical fire behaviour models are based on data collected in experimental fires or in well-documented 
wildfires or prescribed fires (Bachmann, 2001). Empirical models are built by correlating the observed fire 
characteristics with easily measured variables, which describe the so-called fire environment (fuel, 
weather and slope). The process of building up is governed by observations and statistical evaluation of 
forest fires with interrelations taking account of theoretical knowledge (Pastor et al., 2003). 
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Some of the advantages of empirical fire models are: 


The absence of artificiality and scale problems (present in the lab experiments associated to 
physically-based models); 


The integration of numerous factors, which operate in the real world, such as wind and moisture 
profiles (impossible to reproduce in the lab) and fuel heterogeneity. 


The disadvantages of these models are: 


They are basically built upon experimental data and are not to be used in conditions that are 
different to the ones they are meant for; 


Various reasons of economical, practical and other nature virtually limit the possibilities for 
experiments corresponding to real conditions, thus reducing the capabilities of built up models; 


The empirical models thus created are single-dimensional and unalterable. 


According Mendes-Lopes (1998) the future will hopefully bring more fundamental solutions but, for the 
time being, the most apparent benefit of a physical approach to fire modelling is its contribution to 
understand fire propagation mechanisms, therefore helping the experimental design and interpretation of 
field trials. 
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1.4.3 Semi-empirical Models 


Semi-empirical models represent a further development of empirical models and they usually comprise 
one basic equation describing the energy balance and involving the basic (or some of the basic) heat 
transfer methods: radiation, convection and conduction. The heat transfer processes involved are described 
by submodels (Catchpole, 2002), which are largely empirically derived and associated to the type of fuel 
and the region involved (Simeoni, 2002). This is what lessens the universal application of these models. 
Semi-empirical (semi-physical) models are largely single-dimensional or two-dimensional models. When 
such models are to be created it is necessary to establish a balance between universal application 
characteristics and the complex design of the model (Mandel, 2004). 


One of the most common semi-empirical models is the Rothermel model (1992). It is based on physical 
principles but also involves a number of parameters derived from laboratory experiments or from field 
data from Australian grass fires. The physical grounds for this model provides possibilities to use it within 
a wide range of conditions and also to achieve good results in predicting fire behaviour. However, 
according to Catchpole (2002) some of the drawbacks of the Rothermel model include its super sensitivity 
to the height (depth) of the vegetation bed and low efficacy in predicting fire behaviour in vegetation bed 
comprising a large amount of “live” components. 


Another widely used model is the BEHAVE fire prediction system. This is a computer fire simulation 
program that uses as inputs the fuel, weather and topography properties to calculate quantitative fire 
parameters — rate of fire spread, fire line intensity and flame length (Andrews, 1986). Behave equations 
are also based on mathematical and physical laws of thermodynamics and heat transfer. To reflect its 
expanded scope, it is now called the BehavePlus Fire Modeling System (Andrews, 2005). 


Dimitrakopoulos (2001) uses the Behave model to design the Novel nomographs for determination of fire 
behavior characteristics for developed fuel model. Combustion models are made up of vegetation species 
characteristic of the Mediterranean and are considered representative of the region of Greece. 


The following semi-empirical equation are basis for determination of ROS with using of BEHAVE 
software: 


_160+¢, +6) 
PQ 


ROS 


, [m/min] (1.4) 


From this equation can be see that ROS is affected by the next factors: fire reaction intensity I, — rate of 
heat release per unit area of flaming fire front, propagating flux ratio €, fraction of total heat released from 
the flaming front that is absorbed by adjacent unburned fuels and depends on the geometrical properties of 
the fuel bed (surface area/volume ratio of fuel particles, wind factor ®,,, and slope factor ®,, the effects of 
wind speed and slope on the rate of fire propagation, porosity of fuel bed p, effective heating number , 
fraction of total fuel load that is heated to ignition temperature, heat of pre-ignition Q — heat required to 
bring a unit mass of fuel to ignition temperature. The last factor heavily depends of fuel moisture content. 


The rate of heat release per unit length of flaming front (fire line intensity I) is dependent on surface load 
SL and heat content of vegetation fuel H and it is determined using the the following equation: 
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I=SLxHx ROS, [kW/m.s] (1.5) 


The equation used to determine the length of the flame Ly shows that the latter is modelled as a function of 
fire line intensity: 


L, =0.45(1)°”*, [m] (1.6) 


Another commonly known system used to predict the fire behaviour for surface fires is FARSITE (Finney, 
1994). The FARSITE is computer program that uses Huygens’ principle of wave propagation to expand 
fire fronts. In general, Huygens’ principle enables a logical implementation of existing fire behaviour 
models. Each point on the fire front contains information on the time, direction, and rate of fire spread. 
These are essential components of existing models of surface fire spread, fire acceleration, crown fire and 
transition to crown fire, as well as spotting. Fire spreads more rapidly in the direction of the wind and the 
direction of upslope, so an ellipse is often used to quantify the shape of a point source fire. 


FARSITE incorporates models for surface fire spread (Rothermel, 1972; Andrews, 1986) as well as 
transition to crown fire and crown fire spread. The FARSITE model requires the user to identify the data 
files (containing landscape, weather, and wind data). The mouse is then used to input ignitions on the 
displayed landscape. These ignitions can be points or existing fire shapes (drawn as a series of line 
segments). In a similar fashion, users can make minor modifications to the landscape including control 
lines or fuel type changes. The duration of the simulation is determined by either time elapsed or defined 
by the desired ending date and time. The model requires the support of a geographic information system 
(GIS) to manage and provide landscape data. 


The FlamMap fire behavior mapping and analysis system fills a need not met by either FARSITE or 
BehavePlus (Finney, 2006). It uses the spatial fuel and terrain data that is used by FARSITE and 
implements fire models that are in FARSITE and in BehavePlus. The focus of FARSITE is to simulate 

fire growth and the changes that occur over time for a specific fire. FlamMap is a spatial implementation 
of the point models in BehavePlus without the simulation techniques required by FARSITE; each point on 
the landscape is an independent calculation. A map is produced for an area of any modeled value such as 
flame length. Comparisons can be made between locations, or the effect of fuel treatment can be examined. 
Like FARSITE and BehavePlus, FlamMap is a PC based program designed for use by local fire managers. 
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In Catchpole et all. (2002) is presented a semi — physical model for the steady spread of fire through a 
homogeneous fuel bed. The fuel bed is modelled as an arrangement of homogeneous particles all with the 
same moisture content at ambient temperature. The role of the air between the fuel particles is to supply 
oxygen and the thermal capacity of this air is neglected. The diagram of the computational range is 
presented in Figure 1.1. The model is based on a physical representation of the heat transfer processes. 
The heat from the fire flame is transferred to a cell into the fuel bed based on the radiation from the flame 
and the convection resulting from the hot air. The cell in the vegetation layer increases its temperature and 
ignites when it reaches the ignition temperature (a pre-assumed value) and eventually becomes a source of 
heat. The temperature distribution in the fuel bed is derived based on a simple differential equation 
presenting the thermal energy balance for a surface volume element from the fuel bed. When the flame 
reaches the volume element, the temperature reaches ignition temperature and the volume element bursts 
into flame. This is often referred to asan integrated or global energy balance. 


The flame height and flame angle are function of Byram’s fireline intensity (Byram, 1959) and wind speed. 
In cases where the wind velocity in the fire zone is zero, the flame height is equal to its length and this 
characteristic is a function of only the fireline intensity. 


The temperature profile in front of the fire front is modelled by a reducing exponential law of a maximum 
value along the line of the fire front. The maximum temperature is derived based on an exponential 
relation and is a function of the surface area to volume ratio, wind velocity and fuel bed packing ratio. 
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Biomass Properties and Fire Prediction Tools 
This model can be viewed as a bridge between semi — empirical (semi - physical) models and prediction 
tools based entirely on physical laws for the preservation of mass, for the amount of movement and 
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Figure1.1 Schematic diagram indicating position of the fuel element on the surface of the bed. (Adapted 
from Catchpole et all. (2002)) 


Simeoni et all. (2002) present a semi - physical model for determining the surface area occupied by the 

fire (the position of head fire front and rear fire front) as a function of time. The report also presents 
results for a fuel bed comprised of Pinus pinaster needles from the region of Portugal. The basic equation 
presents a non-stationary variation in the vegetation bed temperature and is based on the balance between 


the thermal energy emitted and absorbed by the corresponding computational cell. Heat generation is a 
function of the heat content of the vegetation layer and an exponential law for the variations in the mass of 
the vegetation bed. The parameters of this law have been determined experimentally based on 
intentionally fired fires into a aerodynamic tunnel. The presented model takes into account the effects of 

air flow velocity and the slope of terrain on fire behaviour characteristics. An additional member was 

added to the model in order to take into consideration the effects of the slope of terrain and this member 
represents a radiational heat source and involves a parameter which is a function of the flame tilt angle, 
the emissivity of the flame, the absorptivity of the fuel, and the view factor. This parameter is presented 


using an empirical relation and experimentally determined values for the constants involved in it. 
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Figure 1.2 Diagram of the computational area (on top) and an algorithm for the calculation of the outlines 
of the fire (Simeoni et al., 2002) . 


1.4.4 Physical Models 


Physical models consider various types of heat transfer and recently (in the latest approaches) a multi- 
phase medium is employed to describe the fire. The said medium comprises one fluid phase and a number 
of solid phases and a detailed description of thermodynamic processes and fluid flow (with appropriate 
equations) is provided (Morvan, 2002a: 4-6). The implementation of these models and obtaining accurate 
results requires fast computing equipment, availability of specialised application software (or custom- 
made software) incorporating relevant mathematical models complying with the specifics of the physical 
processes involved, and is highly time consuming (Linn, 2002). 
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The next section of this literature review provides an overview of the major characteristics and 
peculiarities of a fire behaviour model of this type — physical 


1.4.4.1 Modelling of Vegetation Layer 


In some physical models the vegetation layer is presented as a porous media (Zhou et al. 2000, Paz et al., 
1998). In this case the processes considered involve the burning of one solid phase, representing the 
vegetation layer, and one gaseous phase. When no fire is present the gaseous phase is an air medium, and 
when a fire is started into this medium the result is the generation of the products of thermal 
decomposition of the fuel and the subsequent oxidation associated with the burning of the vegetation 
layer. In this particular approach the vegetation fuel bed represents a heterogeneous system made up of a 
solid matrix with a randomly orientated structure. In a small test amount the solid phase coexists with the 
gas phase and a simulation of the processes of interaction between the vegetation particles and the gaseous 
flow is made. The porous medium is in general considered a fluid zone with static solid vegetation 
particles present in it. The porosity of a gas phase is defined as the ratio of the volume occupied by gas to 
the entire volume of vegetation layer. In the similar way the packing ratio of solid phase is the ratio 
between volume occupied by solid phase and entire volume. 


Solid vegetation particles exercise resistance to the gaseous flow and reduce its velocity within the fuel 
bed. The drag forces are described by means of adding a source term to the standard equation for the fluid 
medium. The source term is composed of two parts: a viscous loss term and an inertial loss term. 
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For low-speed flows and laminar flows, the pressure drop is typically proportional to velocity and Darcy's 
Law: 


Vp =—-4% (2.7) 
a 
At high flow velocities, the constant representing inertial resistant factor provides a correction for inertial 
losses in the porous medium. This constant can be viewed as a loss coefficient per unit length along the 


flow direction, thereby allowing the pressure drop to be specified as a function of dynamic head. 


Coefficients accounting for pressure drop are directly dependent on the porosity y. A number of 
dependencies exist to represent the relation between vegetation bed characteristics and these coefficients 
and the major part of them were derived based on experimental data. In Garzon (1997) permeability is 
expressed in the following way: 


Db 3 
a= par ; (2.8) 
180(1—7) 
[D;.f(,).aD, 
D,» = Me (2.9) 
[D;.£(D,).aD, 
0 


In these equations Dp is the average particle diameter and f(Dp) is a function of the distribution L is the 
bed. These dependencies are basically applicable to spherical shape particles (which is different to the 
type of medium being considered here) but they can also be used for a wide range of tasks (Garzon, 1997). 
The following expression is used to describe inertial pressure drop: 


Pini a Cp? DVinggV (2.10) 


In the above equation c; is a friction factor and it is equal to 0.5 (Garzon, 1997). 


An expression for the drag coefficient (Cp) in a bed of pine needles (Pinus pinaster) has been developed 
from Ventura (2002), using as parameters the bed porosity, the airflow velocity and temperature: 


Cy = (—0.1546U cote + 13970/ oso87( 55 + 04913] (2.11) 


This expression is valid for the range of parameters covered by the set of experiments, namely 
Porosity: 0.92 <6 < 0.94 
Flow velocity: 1.15 < U needles [m/s] < 2.48 


Flow temperature: 293 < T [K] < 337 
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Figure 1.3 Pressure drop along the work module, as a function of air 
flow velocity, for 3 different values of porosity (VENTURA et al., 2002) 


This expression can be useful in existing mathematical models for the prediction of surface fire 
propagation. Above results may be used to determine the resistance of the fuel bed. The resulted pressure 
drop data is presented on the Fig. 1.3. 


Morvan et al.(2002) present a model of surface fires in shrubs and grass layer and the vegetation is 
presented via m number of solid phases and one gas phase. Each solid phase is represented as a mixture of 
dry material, water, charcoal and mineral residue (ashes). 


This particular approach of presenting the vegetation layer is also used in the multiphase model described 
by Simeoni et al. (2002). 


1.4.4.2 Mathematical models 


In Morvan et al. (2002) a grouping has been carried out and the basic features of various approaches have 
been presented for the physical modelling of the fire behaviour for fires occurring in a bed of vegetation. 
Major attention has been paid to models based on a multi-phase medium (involving multiple solid phases) 
and also models where the vegetation layer is considered as single-phase discrete medium (and in some 
cases a porous media). 


In the approach presented by LARINI et al. (1998) PORTERIE et al (2000), MORVAN et al (2000, 2001, 
2002, 2003) the combustible medium is considered as a multiphase medium composed of a gaseous phase 
and several solid phases. Each solid phase consists of the same fuel particles, which have the same shape, 
the same size, the same physico-chemical properties, thus the same behaviour with respect to a fire. In 
particular, each family of particle is characterised by its surface-to-volume ratio and in order to achieve a 
spatial description of the medium, the fraction of space volume occupied by each family or each solid 
phase, at each point and each time, must be known. 


In the approach of SERO-GUILLAUME et al. (2002), the vegetal phase is considered as a fractal medium, 
which can be described mathematically. Moreover the different parts of plants are considered as a porous 
medium (Zhou et al., 2000), in order to take into account the internal structure of vegetation. 
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When the fire approaches the particle it receives energy from this fire by convection and radiation. Its 
temperature is then raised from the ambient temperature up to the boiling water temperature, at which it 
looses its water. Subsequently, the temperature rises again and the plant material is decomposed step by 
step releasing pyrolysis gases. The reactive part of these gases is combined with the oxygen in the air 
(gaseous phase chemical reactions). The gas temperature then increases, causing its expansion and, due to 
buoyancy forces, the gas moves. These movements play a very important role in the necessary transport of 
oxygen and also in the energy transfers. At the end of the pyrolysis, the particle is mainly composed of 
chars. Oxygen may reach the surface of the particle and react with the remaining chars. It’s the 
combustion of chars, which causes the regression of the particle surface. 


The description of the medium is based on determining the parameters of the gaseous phase velocity, 
temperature, mass fraction of chemical species and pressure, as well as on determining temperature, mass 
fraction of chemical species and density at each point of the motionless particle. 
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At this scale, in order to calculate the state of the medium at each point of the gas, the balance equations of 
reactive fluid mechanics are available (mass, momentum, energy and chemical species) and at each point 
of the particle, a system of balance equations mass, chemical species, energy) is available, which describes 
the behaviour of a solid conductive medium that can loose its water and be decomposed into pyrolysis 
products. In addition, at the interface between the solid particle and the gas, some conditions must be 
verified (local interface conditions). 


Solving all these equations at this scale can be done only in very simple geometric configurations, but 
cannot be done in real configurations with numerous particles distributed in a gas. In this last case, and at 
this scale, it is impossible to obtain solutions of the strongly coupled balance equations at each point of the 
gas, each point of each particle, which in addition must verify all local interface conditions. 


Thus, it’s necessary to change the scale of observation, in fact, to increase this scale. From a mathematical 
point of view, operating this change of scale is equivalent to average the equations established at a point of 
the gas or at a point of the particle over some volume around these points. It’s a usual operation for 
modelling multi phase multi constituents materials and several methods of averaging are available. 


After some mathematical transformations, a system of equations is obtained for each phase (gaseous phase 
and solid phases). These equations establish the relations in time and space between the weighted average 
of the different variables (temperature, velocity, density,...). It is worth noting that the averaging operation 
leads to the appearance of new terms in the balance equations, which are due to the interactions between 
the different phases (here between the gaseous phase and each solid phase). Thus, systems of partial 
derivative equations are obtained that are strongly coupled through these terms of interaction. For instance, 
these terms of interaction are due to the mass transfers related to the drying and pyrolysing process of the 
particles, to the drag forces, to the heat transfers by convection and by radiation. system of multiphase, 
reactive and radiative equations, rigorously deduced from the instantaneous point equations (generalised 
NAVIER-STOKES equations for the gas and point balance equations for the solid), is now available. 


The interaction between the continuous and the discrete phase is described mathematically in several basic 
ways. Source terms are added to the equation describing the preservation of the mass of the gaseous phase 
and these are used for the modelling of the moisture evaporation (in the form of water vapour) and the 
release of volatiles from the vegetation layer. To maintain the mass balance for the entire area being 
considered, a variable mass of the solid phase is also modelled. Moreover, some authors also assume 
constant density of the vegetation particles throughout the pyrolisis process (Zhou et al. 2000), while other 
models consider a reduction in the density of the solid phase as a result of moisture evaporation and 
release of volatiles, and the volume occupied by this solid phase stays constant (Marcelli et al., 2004). 
Turbulence processes in the continuous medium are modelled mainly based on the k — ¢ turbulence model 
or RANS models, but there are also other approaches, such as using large eddy simulation (LES) 
(McDonough et. al., 2004, Hostikka et al., 2001). 
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1.4.4.3 Modelling heat transfer 


With the physical approach used for modelling fire behaviour it is especially important to consider the 
various processes of heat transfer within the gaseous medium and the vegetation medium both 
individually and between them. These types of models also give consideration to the three types of heat 
transfer (conduction, convection and radiation) and it is especially important for the spreading of the fire 
to consider the transfer of hot air from the fire zone and the heat radiation from the fire flame. The heat 
transfer via conduction heat transfer is of significance for vegetation layers of high bulk density (forest 
litter) in the direction of the depth of the vegetation layer. 


Zhou et al. (2000) express the thermal interaction between the gaseous phase and the solid phase for pine 
litter using the following convective/conductive heat flux: 


Q, =Ah,(T-T,) (1.12) 
In the above relation A is specific wetted area and this is proportionate to the surface to volume ratio of 


the solid phase, T u T, are the temperatures of the gaseous phase and the solid phase, respectively. The 
heat transfer coefficient is deduced from Nusselt number of solid phase 


Nu = ne = 0.683 Re’** (1.13) 


Ss 


In above equation A, is thermal conductivity of solid phase. 


The important mechanism for heat transfer from the flame zone is radiation. The radiation heat flow from 
the high temperature zones is a function of the absorption coefficient of the medium a, the Stefan — 
Boltzmann constant o, and is proportionate to the fourth degree of the temperature of the emitting body. 
For the particles inside the vegetation layer located close to the fire flame the increase in temperature is 
mostly dependent on the amount of heat they absorb as a result of the radiation heat transfer. The 
expressions for the radiation heat sources for the gaseous phase and the solid phase are determined as 
follows (Zhou et al., 2000; Johnson et al., 2001): 


04, = a (Aor =G) (1.14) 
Qras = §.a,(G —40T:) (1.15) 


The indices g and s in (1.14) and (1.15) above are used to indicate the parameters of the gaseous phase and 
the solid phase. The absorption coefficient for the vegetation layer is calculated using the standard formula 


a, =A/4 (Porterie et al., 1998). For pine litter featuring high values for wetted area A the thermal flow from 


the flame radiation will affect the temperature rise inside small depths of the vegetation layer. Thus, a 
conclusion can be drawn that the convection heat transfer and the radiation heat transfer from the fire zone are 
especially important for the increase in the temperature at the surface of the vegetation bed. Major role for the 
temperature rise of the vegetation layer at locations close to the fire flame plays the heat radiation, and as the 
distance from the fire flame is increased, the effect of the hot gas convection becomes more significant. 
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The absorption coefficient of the gas phase is related to the concentration of combustion products (CO2, 
H,0) and soot formation (Kaplan et al., 1996). According to Marcelli et all. (2004) “The chemical nature 
of the burning material and local specific fire conditions can lead to a large production of soot” and also 
“An accurate prediction of radiative transfer requires accounting for soot production/destruction, 
especially in turbulent flames”. 


1.4.4.4 Modelling the heat generation 


For the burning process to occur and persist, it is necessary to have heat generation in order to heat up the 
surrounding vegetation and cause the release of combustible gases, which release additional amount of 
heat as they undergo further oxidation. Thus, two steps are identified during the process of heating up the 
vegetation layer: the first one is a endothermal process of heating up and evaporation of moisture and 
release of volatiles, leading to subsequent oxidation accompanied by exothermal heat generation and 
continuing temperature rise. 
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In the course of the modelling two groups of models can be identified depending on the presentation of the 
heat release in the process of burning of the fuel. In the first group heat is generated once an ignition 
temperature (experimentally determined) is reached in the corresponding cell in the vegetation layer and 
this heat forms a heat source included in the energy equation (McDonough et. al., 2000; McDonough et. 
al., 2004; Paz, 2000). This heat source is presented as a function of the heat content of the vegetation 
particles and the rate of variation of the vegetation mass. Normally, the ignition temperature is assumed to 
be within 300 - 400°C (after the moisture in the fuel has been released) and therefore, the value assumed 
for the amount of heat to be released during the burning process is the value which corresponds to the 
higher heating value (HHV). The variations in the mass of the fuel result from the thermal decomposition 
which is in turn dependent on the rate of heating up of the vegetation, as well as on the temperature and 
the chemical composition of the vegetation layer. Very often variations in the mass are modelled based on 
the rate of devolatilization of the fuel. Due to the lack of sufficient data for the different types of 
vegetation involved, the kinetin rate is used with values for the activation energy and pre-exponential 
factors based on data for coal particles (Paz, 2000; Zhou et al., 2000; De Bruyn Kops et al., 2004). 


The second group of models considers the processes of thermal decomposition of the vegetation and the 
oxidation of volatiles, char and tar. Modelling of chemical processes is used and these involve one to 
several chemical reactions (Zhou et al., 2000; De Bruyn Kops et al., 2004; Marcelli et al., 2004) and they 
only provide one very general presentation of the complex phenomena involved in the combustion of the 
vegetation layer. The right choice of the system of chemical reactions, the rate and the type of the 
chemical processes, as well as the rate of release of the reagents and their chemical composition, are all 
highly significant for the amount of heat released in the fire zone and hence, for obtaining the physical 
picture of the temperature distribution, of the distribution of the other parameters of the gas flow, as well 
as the characteristics of the fire behaviour. 


1.4.4.5 Modelling the mass loss rate of vegetation layer 


The rate of conversion of the vegetation into gaseous products maintaining the burning process is not 
among the characteristics of fire behaviour but it affects to a considerable degree the intensity and rate of 
spread of the fire. It is a well known fact that reductions in mass of solid fuel are associated with the 
heating up and the accompanying processes of moisture evaporation, release of volatiles, and formation of 
tar and solid carbon. The latter three undergo oxidation in the course of the burning process. It is exactly 
the rate of generation of these products that is expressed by means of mass loss rate. The main factors that 
affect the reduction of the mass of the vegetation bed in the course of thermal decomposition are the rate 
of temperature rise and the fuel temperature. The method used for determining mass loss rate is based on 
the differential thermal analysis (DTA). However, when this method is used the rate of temperature rise is 
rather limited and therefore, it is possible to obtain more reliable data for the particular characteristic being 
considered here only when specifically designed laboratory devices are used to study the fire behaviour 
for intentionally initiated fires. 
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Figure 1.4. Time evolution of the rate of mass loss (Fuel + PP, D = 0.4 mm, Plate = Yes) 
(Dupuy et al., 2002) 


In their reports Leoni et al., (2002) and Dupuy et al., (2002) describe experimental rigs and present 
measurement results on mass loss rate for pine litter. These experiments were conducted with the purpose 
of creating models for fuel mass variations throughout the burning process. This data can be useful in the 
modelling process: fuel moisture evaporation rate, combustion gas emission rate, and also in comparing 
these results with the results of the numerical modelling of a process of fire spreading in a vegetation bed. 


1.5 Thermal Decomposition and Combustion of Vegetation 


1.5.1 Thermal decomposition of vegetation 


What are the changes that occur in vegetation when heat is transferred? In short-term heating of up to 
100°C no chemical changes are observed in wood. The moisture content of wood is evaporated. Heating 
wood up to 125°C completes the process of evaporation of the moisture content. When wood temperature 
reaches up to 210°C, chemical reactions take place in wood as a result of which a mixture of gases is 
generated (70% carbon and 30% carbon-based combustion mixtures), along with some volatile 
compounds (CO, H etc.), thus reducing the mass of wood. Ifa flame is brought near wood at this 
temperature it bursts into short flames, i.e. the wood ignites. This temperature is called ignition 
temperature. Heat generation occurs as a result of the process of oxidation of wood and the reaction is 
exothermal. At temperatures of 260°C wood combustion is steady-state and the temperature is referred to 
as flame-burning temperature. As the temperature of wood increases the amount of combustion gases 
generated increases and the generation of non-combustion CO, reduces. At temperatures of 450°C the 
process of generation of gases from wood stops, flames disappear and surface combustion of the coke 
residue (carbon) is initiated, which rises the temperature up to 900°C (Videlov, 1993) 


Download free ebooks at bookboon.com 


39 


Please click the advert 


Biomass Properties and Prediction Tools for 
Biomass Properties and Fire Prediction Tools Vegetation Wildland-Urban Interface (WUD) Fires 


Fuel oxidation can be described from the following simplified reaction: 


Wood + air => water + carbon dioxide 
(CoH iOs)n +6n O, => 5n H,0 +6n CO, 


The stages of vegetation oxidation are (Johnson et al. 2001): 


e Ignition — wood oxidation occurs in a single point of the wood surface; 
e Inflammation — oxidation spreads along the entire surface of wood; 
e Combustion — oxidation normally spreads on the external surface of wood. 


Wood undergoes the following stages throughout the burning process — evaporation of moisture content, 
evaporation of volatile compounds, generation of fluid products, combustion of hard coke residue and 
glare (Byram, 1959a: 64-69). 


Wood consists approximately of 50% carbon, 44% oxygen and 5% hydrogen | by weight (Johnson et al., 
2001: 57). 5,9 kg (4,6 m°) of air is required for the complete burning of 1 kg air-dried wood. The mineral 
residue of the combustion of wood is 0.5% of its initial mass (Videlov, 1993). 
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The Drying is energy consuming (exothermic) process and is considered to be finished when the moisture 
contents of the fuel are evaporated. The moisture contents, achieved by proximate analysis are definitive 
for the description of the entire process. The heat-up of the fuel is evaluated as the rate of the change of its 
temperature, i.e. [K/s]. At low drying rates (usually 10 — 10° K/s) the water contents are diffusing to the 
fuel surface and evaporate, as well as released through the pores of the sample. It worth saying that the 
structure (ash lattice) of the fuel sample is definitive for this process. At these levels of heat exchange the 
temperature of the sample could be considered constant and homogeneous. At high heating rates (up to 
10° — 10° K/s) the moisture contents are released rather intensive and even could evaporate inside the 
pores of the fuel particle (no homogeneous temperature distribution in the fuel particle volume is 
observed), thus increasing internal pressure and even burst the solid matter. However this rarely happens, 
so one could assume more convenient (especially for modelling purposes) process, which occurs at low 
heating rates and resulting homogeneous temperature distribution. At drying regime the temperature of the 
fuel sample is relatively low (up to 100°C at atmospheric pressure) and no intensive thermo-destruction of 
the hydrocarbons proceeded. The energy, consumed for water evaporation in case of forest fires is energy 
(predominantly supplied by radiation) released during combustion of the biomass fuel and rules as 
extinction factor for the burning process. 


Devolatilization. When the drying process has finished the biomass fuel consists only of combustible 
matter and ash (mineral matter). The combustible matter consists of volatiles, tar and solid carbon residue. 
The volatile gases are products of thermal degradation of the combustible matter. Due to external heat 
supplied to the fuel the chemical chains of the complex hydrocarbons are broken and light organic matter 
— volatiles, is produced as well as tar. There is deep commitment between the heat-up rate of the fuel 
sample and yield of volatiles — as the heating rate is increased the released volatiles quantity is higher and 
respectively the tar yield is relatively lower. The thermal degradation is basically exothermic process and 
three staged of this process could be observed (Johnson et al., 2001). The energy, consumed for light gases 
devolatilization is called latent energy of devolatilization. Some authors (Kjaldman, 2000) ignore this 
energy, when come to build model for the devolatilization process, but principally this heat quantity is 
“hidden” .into the latent heat for water content evaporation. 


Experimental observations as well as researchers data show that at thermal degradation (pyrolysis process) 
around 80% of the mass is released through volatile gases (Richard et al., 2002). This again confirms that the 
vegetation material could be treated as gas releasing surface in close mass and energy interaction with the 
homogeneous phase. The remaining solid phase release tar and char residue. The released tar is secondary 
thermal degradation based step and is converted into gas phase by oxidation. The char residue is consumed 
by surface reaction, which is diffusion limited, i.e. the diffusion of the oxygen to the reacting surface is the 
limiting factor and thus the reaction determining one. The experimental investigations of the combustion 
process of wood based biomass show that the volatile gases are released in the form of non-oxidized 
chemical species in the gaseous phase, primarily as CH4, CO, H2, O2, C,H, and these components undergo 
chemical reaction in the following mixing process if appropriate conditions occur (Morf et al., 2002). 
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The pyrolysis temperatures for the three main constituents of the combustible matter of vegetation vary as 
follows: hemicelluloses at 200-260 °C, cellulose at 240-350 °C and lignin at 280-500 °C. The extended 
period of contact between the primary products allows their re-polymerization, cracking and/or cross- 
linking to produce char. Consequently at high heating rates the yields of char decreases (20 % of original 
cellulose weight at about 300 [°C/s], compared to about 2% at 1500 [°C/s]). With an appropriate heating 
rate thin samples of cellulose (as pine litter can be considered) can be completely volatilized with no 
production of residual char. In principal besides of the heating rate, char formation is a function of particle 
size and the moisture contents as well as the mineral matter properties and ash contents. The experiments 
show that with increasing amount of inorganic contents the char yield is increasing. Thus for relatively big 
fuel particles (or agglomerations of the pine litter matter) and increased moisture content and at slow 
heating rates the result is high char formation and relatively low combustible matter conversion. 


The three-step pyrolysis process is described in the following diagram. 
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Figure 1.5 Three-step reaction scheme for the pyrolysis process 
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The consumption rates are Arrhenius type expressions as follows: 


production rate of the volatile gases k, 


k, =3.17x10" exp( A [1/s]; (1.17) 
Uu. 
production rate of the char + tar k, 
k= 1.3200" ea x [1/s]; (1.18) 
RuT 


In equations (1.17) and (1.18) R, is universal gas constant (Ru=83 14 [J/(kmol.K)]) and T is the absolute 
gas temperature, [K]. 
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The presented values for the activation energy and the pre-exponential factor given here are corrected and 
adjusted as the values are specific for each reaction and are function of the environment conditions of the 
process (Peters et al., 2002; Richard et al., 2002; Brush et al., 2003). 


1.5.2 Combustion Process 
1.5.2.1 Turbulent combustion models 


Principally the problem of the mathematical modelling of turbulent-combustion processes is reduced to the 
task of modelling the interaction between the turbulence and the chemical reactions. The combustion process 
in forest fires occur at conditions of continuous changes in the temperature, concentration of chemical 
components, velocity and chemical state of the reacting species. Under these conditions the process can be 
presented mathematically by formulating a complex system of non-linear differential equations, which 
covers following areas: flow approximation, turbulence, chemistry, heat transfer (radiation, convection and 
conduction), thermochemical data for the species involved, problem geometry and initial and boundary 
conditions, and application of complex mathematical model to involve all these models. 
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The determination of the chemical reaction rate is the basic problem in combustion processes modelling. 
There is a large variety of models developed for turbulent combustion modelling. There are more then 100 
models known and published in the literature. Their description is a huge task, but the models most widely 
used for practical predictions can be grouped in several approaches and classified according their 
development, application and improvement. Some of these models will only be mentioned in the following 
section and the emphasis will be on applied combustion model — eddy dissipation model (EDM). 


Turbulent combustion models, based on the “flamelet theory”. These methods assume that the turbulent 
reacting flows consist of an ensemble of reaction-diffusion layers that are continuously displaced and 
stretched within the turbulent medium. The structure of the flamelets can be approximated by employing the 
flame sheet model or by assuming equilibrium or partial equilibrium. Several “flamelet” based models for 
turbulent combustion are developed. Some of them are: Flamelet Surface Density Models, Algebraic model 
of Bray-Libby- Moss, Model using transport equation for the mean flame surface, etc. 


Probability density function models (PDF). Taking into account the stochastic character of the turbulence it 
is logical to use the Probability Theory for mathematical describing of the turbulent flow field. Further the 
method is extended in deriving a transport equation for a joint (for the velocity components and for the 
scalars) probability density function. Thus becomes possible to define the velocity and the scalar 
components together. Therefore the joint PDF (JPDF) approach is widely used in turbulent combustion 
modelling. The idea is that the realistic chemistry, which includes a detailed chemistry, can be presented 
clearly, without any assumptions for the flame structure. However, the effects from the molecular processes 
and the pressure fluctuations gradient appear in the equation as conditional expectation, which has to be 
modelled. Obviously, the transport equation includes many unknown terms. The solution of the JPDF 
equation increases drastically the CPU time, thus leading to problems with the numerical implementation. 


The turbulent models based on the eddy dissipation approachs are widely used in combustion modelling 
society because of their capacity to represent the combustion physics, computationally not so expensive as 
well as achieved results are qualitatively comparable and close to the experimental data. 


Eddy Break Up Model (EBU). This model was presented by Spalding and widely used in the 1970s. The 
idea of the EBU model is based on the following assumptions: combustion proceeds as an infinitely fast 
reaction; the chemical reaction rate R, is proportional to the rate of dissipation (¢/i) of the individual 
eddies; the whole process is presented as a single, irreversible one-step reaction: 


[1kg of fuel] + [Veue kg of oxidizer] — [(Vrue+ 1)kg of products] (1.19) 


These assumptions lead to the following mathematical expression for the reaction rate Regu : 


é ! 
Regu = Cray P— - - , [kg/(m’s)] (1.20) 
where Cggu is model constant, Vpy1- stoichiometric coefficient (amount of oxygen ,required to oxidize lkg 


of fuel at stoichiometric condition ), « - turbulent kinetic energy , [m’/s’], ¢ - dissipation rate of the 
turbulent kinetic energy , [m’/s*], p — gas density , [kg/m’]. 
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This equation expresses ideally the fact that the local chemical reaction rate depends on the gas mixture 
fluctuating components Y°. However, in practice this leads to a number of complications related to 
defining almost unpredictable variables such as fluctuations. The application of the EBU model is also 
mostly limited to modelling only pre-mixed combustion [PhD thesis Dobrin]. 


1.5.2.2 Eddy Dissipation Combustion Model (EDM) 


This model is proposed by Magnussen and Hejertager (1977). The model develops the conception of eddy 
dissipation and overcomes successfully some problems related to EBU. In this model the chemical 
reaction rate is presented by the mean mass fraction values of the reacting species instead of their 
fluctuations — actually this is the main advantage of EDM compared to EBU. The main factor influencing 
the rate of infinitely fast chemical reactions is the diffusion of species. The chemical reaction itself takes 
place when the reagents have mixed up to the molecular level at high temperature. The chemical reaction 
rate is determined by the mixing rate, i.e. by the rate of dissipation of eddies. The fuel and oxidant appear 
as fluctuating intermittent quantities so that a correlation is possible to exist between the fluctuations of 


: : Yuvy" F 
the chemical components and their mean values (7 ~ ). Thus the reaction rate REDM can be 
expressed by the eddy decay rate and average mass fraction of the reacting elements. 


Because of its simplicity and testified applicability it is one of the most popular concepts for modelling 
local reaction rates in practical turbulent combustion systems. This model has the ability to handle 
premixed, partially premixed and non-premixed combustion with sufficient accuracy. 


1.5.2.3 Methane oxidation model 


Methane oxidation will be reviewed in order to present developed models for chemical reactions of CH4 
as representative for hydrocarbons. This oxidation model is used as basis for development set of chemical 
reaction that is applicable to complex problems, such as biomass and vegetation combustion process. 

The process of methane oxidation is well known and therefore several mechanisms for the kinetics can be 
found in the literature. These mechanisms can be classified as: 


e Full mechanism describes the interaction between hydrocarbons and the so called nitrogen 
chemistry. It consists of 225 elementary reactions and 48 participating species. Detailed 
description of each of these reactions is not necessary when practical engineering problems are 
considered, considering that their numerical investigation is limited by the computer possibilities. 

e Skeletal mechanism is obtained from the full mechanism after analysing the rate of generation and 
consumption and the dependence of the solution to the boundary conditions. Actually the skeletal 
mechanism represents a compromise and balance between complexity, accuracy and applicability 
of the model. 

e The reduced mechanism represents the methane oxidation under ideally stirred reactor and 
therefore it is of special interest for the current investigation. 


The chemical reaction included in four steps reduced mechanism are: 


CH, +2H+H,0 + CO+4H, (1.21a) 
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CO + H,O — CO, + Hp (1.21b) 
2H+M—>H+H+M_ (1.21c) 
0, + 3H) > 2H + 4H,O(1.21d) 


This mechanism is more detailed than one and two step reaction mechanism, but actually in some cases 
makes the numerical modelling unstable and slow converging, the results are close to the less complicated 
and more computationally stable two step reaction model without paying the efforts for describing the 
complex chemistry phenomena. 


Two steps reduced mechanism is proposed [Dryer et al. | and it is valid for high temperature oxidation of 
CH, and CO, as follows : 


ne 
CH, + 3/20) > CO +2H,0 (1.22a) 
Re 


CO +1/20, > CO, (1.22b) 
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The numerical simulations based on the two steps mechanism, depicted with equation (1.22a) and (2.22b), 
show a good approximation of the combustion process. This mechanism can be used for both: initial 
calculation followed by more detailed calculation of the turbulent combustion and for the complete 
calculation of the combustion process. Furthermore, the results are quite close to the experimental data 
and thus reliable applicable in most of the engineering combustion problems. 


One step global reaction mechanism. This model is based on the assumption for an infinitely fast reaction 
where CH, converts directly to CO2 and H,O vapour without considering intermediate species. 


CH, + 20> — CO, +2H,O (1.23) 


The use of this mechanism leads to several advantages when engineering problems are considered and 
therefore it is widely used. However when local flame characteristics are of importance the model does 
not provide the adequate accuracy. Actually it commonly over predicts the temperature field and thus 
introduces significant error in the numerical investigation. 


1.5.2.4 Modelling of the turbulence-chemical reactions interaction 


The inclusion of the kinetics into the calculation is related to the choice of appropriate model for the 
turbulent combustion. There are several known techniques for coupling the turbulence and the chemical 
kinetics, based on the Eddy Dissipation Approach. 


The smallest rate is the limiting. When the turbulent mixing and the chemical kinetics are considered to 
happen simultaneously, it is considered that the limiting factor will be the smallest rate, i.e. the time for 
the reaction is determined by the time needed of the slowest of turbulent mixing or chemical kinetics 
phenomena. Thus the effective rate will have the following expression 


R egective MIN(R poutence > hen eiaies) ? [kg/(m’s)] (1.24) 


where Reftectivee 18 effective rate of the chemical reaction in turbulent reacting flow. 


The series process approach is appropriate for first order kinetic rates. The reacting species should mix 
due to turbulence and once they are well mixed the kinetics are considered. To eliminate the initial 
composition in the well mixed zone, the two processes, mixing and reaction are considered in series and a 
global rate is defined as follows: 


1 1 1 4 
= + , [kg/(m’s)] (1.25) 


effective Taine chem.kinetics 


This approach is testified and gives more realistic results and is widely applied in applications involving 
combustion modelling. That’s why series approach will be used in development the model for forest fires. 
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1.6 Experimental Data for Fire Behavior Modelling Validation 


The fire progression is pointed out as a determinant parameter when evaluating model performance 
(Miranda et al., 2005; Dupuy, 1997; Valette, 2000). This process is described by the main fire behaviour 
characteristics — fire front position, rate of fire spread, flame characteristics (length, height and angle). 
Temperature and velocity distribution are also included as a characteristics of fire progression. 


The second group of experimental data that can be used as a determinant parameter for model validation is 
concentration of pollutants in a given spot also called air quality data. The experimental results from this 
measurements are utilized by smoke dispersion models as a way of validating the performed simulations 
(Miranda et al., 2005). 


Experimental data presented in this section are mainly derived from laboratory experiments but also from 
observations over intentionally caused and controlled fires in pine needles litter. The results presented 
underneath are experimental data for ROS and flame characteristics, and also data for temperature and 
velocity distribution in and above the vegetation layer for vegetation types occurring in the regions of the 
Mediterranean countries: Portugal, Spain, France and Italy. This information is used to analyse the 
influence of various factors and also for the assessment of the variations of quantitative characteristics 
determining fire behaviour, as well as in the process of preparation and planning of experimental work 
involved in determining these characteristics. It is also useful for the theoretical modelling of quality 
based on the formation of various submodels and also as validation data for the results derived from 
numerical experiments. No data for experimental fires in a bed of vegetation originating in Greece was 
found in the course of the research study carried out in available publications on this issue. 


The presented results are mainly for a vegetation bed comprising Pinus pinaster pine needles. 


1.6.1 Rate of Fire Spread and Flame Charakteristics 


Simeoni (Corsica — France 2002) describes an experiment conducted in a wind tunnel under various wind 


velocities (0 — 3 m/s), slopes of terrain (up to 30°) and 10% fuel bed moisture content. Derived data is 
presented with no provision made for considering the values for fire behaviour characteristics, and only 
taking into account the location of the fire front in a given moment after the start of the fire. In the 
particular experiment the temperature was measured at 3 cm above the fuel bed. The maximum measured 
temperature (as indicated on the plots) was about 950°C. 


According to data presented by Zhou (Portugal 2000), for wind velocities of 0.8 m/s ROS = 0.63 m/min, 
and for wind velocity of 1.8 m/s ROS = 0.84 m/min. The presented experimental results were obtained 
from the burning of pine needles Pinus canariensis in crosswind conditions. The vegetation fuel bed was 5 
cm high. The characteristics of Pinus canariensis are: density- 400 kg/m’; fuel load - 0.55 kg/m’; surface 
area to volume ratio - 5600 m’'; cellulosic volatiles - 64.8% and initial moisture content - 0%. 
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Mendes-Lopes (Portugal 1998) reports on the results obtained in several sets of experiments, which where 
put forward at different scales, from field to small laboratory rig experiments. A bed of Pinus pinaster 
needles is burned under controlled conditions (wind speed and direction, terrain slope, and fuel moisture 
content). ROS and geometrical characteristics of the flame are obtained from video records. Results show 
that ROS increases steeply with wind speed for wind-driven fires but does not depend on wind speed for 
back-wind propagation. ROS increases slightly with slope for up-hill propagation, and is not slope 
dependent for down-hill cases. ROS decreases when fuel moisture content increases. Flame angle and 
flame height are also dependent on wind velocity, terrain slope, and fuel moisture content (Mendes-Lopes 
et al., 1998). 


The fuel bed consists of a layer of Pinus pinaster needles. It tries to reproduce a typical L layer found in 
Portuguese stands, with a load of approximately 0.500 kg/m’ on dry basis. The conditioned pine needles 
are scattered uniformly in the tray with a load of approximately 0.500 kg/m” on dry basis, and a fuel bed 
of roughly 3 to 4 cm depth. 


The parameters of the experiments are (Mendes-Lopes et al., 1998): 


e Wind velocity: -3, -2, -1, 0, 1, 2, 3 m/s; 

e = Terrain slope: 0°, 5°, 10°, 15°; 

e Slope orientation: up-slope, down-slope; 
e Fuel moisture content: (1041)%, (18+1)%. 


life- 
changing 
careers 


Where your talents benefit you, and millions of others 


At Novo Nordisk, we are working toward a future where everyone's potential can 
be fulfilled. What about yours? Right now, we are looking for outstanding candi- 
dates to join our 7 Graduate Programmes. This is your opportunity to dream big, 
realise your potential and help defeat diabetes worldwide. 


7 International Graduate Programmes 

Graduate Programmes are available within Global Finance, Global Marketing, 
Supply Chain, Business IT, Research & Development, Business Processes and 
Pharma Management Education. Most programmes are based in Denmark and 
Switzerland with rotations abroad. All programmes start 1 September, 2009. 


Read more about our Graduate Programmes at novonordisk.com/graduates. 
Deadline for application is 12 February, 2009. 


www.novonordisk.com/graduates @ 


novo nordisk” 


Download free ebooks at bookboon.com 


49 


Biomass Properties and Prediction Tools for 
Biomass Properties and Fire Prediction Tools Vegetation Wildland-Urban Interface (WUI) Fires 


ROS increases with the wind velocity at an increasing rate. The effect of slope further enhances this increase. 
Values range from 0.20-0.25 cm/s for no-wind propagation up to 4.7-7.3 cm/s for 3 m/s wind. For no wind 
situation, flame height varies between 25 and 35 cm with no apparent correlation with terrain slope. For 
back-wind flames a small increase of wind speed from zero originates a deep decrease in flame height to 
about 10 cm; further increase in wind speed does not cause any changes to flame height. A different 
behaviour can be observed for wind-driven flames. Flame height steeply increases initially, reaches a 
maximum of about 45 cm for wind speed of, say, 1.5 m/s, and then decreases (Mendes-Lopes et al., 1998). 


Presented in Mendes-Lopes (2003) are some additional results to the study published in Mendes-Lopes 
(1998) and some statistical data processing has been carried out using an analysis of variance (ANOVA). 


Viegas (1998) present results from experiments involving the determination of rate of spread and flame 
characteristics (height and length) for various values of surface load, fuel bed depth, fuel bed porosity and 
fuel moisture for Pinus pinaster pine needles. The wind velocity and slope of experimental table surface 
are of zero value. The results illustrate the fact that for vegetation bed moisture content values higher than 
65% the rate of spread of the fire is equal to zero - it is not possible to maintain the fire in such conditions 
and it goes out. As moisture content is reduced the ROS increases and this increase is significant in a fire 
occurring in a vegetation bed with less than 25% moisture content. According to this publication source 
the maximum values for ROS are 0.4 — 0.5 cm/s for moisture content of 5 %. The fitted curves to 
presented experimental data are third order polynomials. The relation between ROS and surface load is 
linear and the rate of spread of the fire increases as the mass of the vegetation bed is increased. The 
presented results are approximately between 0.1 cm/s for 0.5 kg/m° and 0.25 cm/s for 1.5 kg/m’. The 
length of the flame is presented as a function of recorded values for ROS (0.1 — 0.4 cm/s) and it changes in 
a linear manner within a range of between 2 and 7-8 times the depth of the vegetation bed. 


The literature survey carried out shows that the basic conditions which influence fire behaviour to the 
largest extent are wind velocity, slope of terrain and fuel moisture content. It is possible to provide a 
quantitative assessment, presented analytically, of the influence of these parameters over fire behaviour 
characteristics based on a statistical model developed using a regression analysis of experimental data and 
selecting a suitable analytical presentation for a model equation. In his study Fernandes (Portugal 2002) 
establish quantitative relationships between important fire behaviour descriptors and its environment that 
are applicable to low-to-moderately intense fires in maritime pine (Pinus pinaster) stands. The 
experimental burning program was conducted in six study sites within an elevation range of 450 - 970 m. 
From the statistical analysis of results Fernandes (2002) derived the following equation for ROS 
applicable to pine litter: 


ROS =0.603U°*® exp(—0.035M ,, + 0.0585), (1.26) 


where U — wind velocity, km/h; 
Mgg — moisture content of the surface dead fuel, %; 
S— slope, %. 


Equation (2.26) provides a model for the effect of wind using the power law function, and the effects of 
the slope of terrain and vegetation bed moisture content are modelled based on an exponential function. 
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It should be noted that, although the general behaviour observed is consistent and somehow repeatable, 
there is an appreciable amount of data scattering that makes it difficult for point-to-point or even case-to- 
case comparison with numerical simulations. 


Fig. 1.6 illustrates a comparative plot presenting the results obtained for ROS for pine needles fuel bed 
provided by abovementioned sources. 
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Figure 1.6 Data for ROS variation by various conditions — wind velocity, slope of the 
terrain and moisture content 


Figures 1.7 through 1.13 provide a graphical illustration of experimental results obtained by Mendes- 
Lopes (1998), whose report provides the largest amount of data on fire behaviour characteristics. 
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Figure 1.7 Influence of wind velocity on rate of spread, fuel moisture content (fmc) 10 % 
(Mendes-Lopes, 1998) 
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Figure 1.8 Importance of the fuel moisture content on the influence of wind velocity on the rate of spread 
(Mendes-Lopes, 1998) 
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Figure 1.9 Typical examples of results with different scattering on the discrete values of flame angle 
measured for two experimental runs: a) low scattering (0 m/s, 0°, 10 %); b) high scattering (1 m/s, 0°, 
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Figure 1.10 Influence of the wind velocity on flame angle, fmc 10 % (Mendes-Lopes, 1998) 


50 > 
[ e 

E 40 + ¢?¢ 
= [ o¢ 
~ 30+ o ¢ o 
= a x 7 2 Senne . 
Z o+ e o¢ ~~, 
® | ® 
E o 
z 10 + 

0 


0 20 40 60 80 100 120 140 160 
Time (s) 


Figure 1.11 Typical example of variation of flame height with time, wind velocity 3 m/s, slope 0°, fmc 18% 
(Mendes-Lopes, 1998) 
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Figure 1.12 Influence of the wind velocity on flame height, fmc 10 % (Mendes-Lopes, 1998) 
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Figure 1.13 Importance of the fuel moisture content on the influence of wind velocity on flame height 
(Mendes-Lopes, 1998) 


1.6.2 Temperature and velocity distribution for and above pine needles litter 


Reports provided by Ventura (1998), Marcelli (2002 and 2004) present measurement results for 
temperature in and above the fuel bed. The maximum temperature is within the range of 1000 — 1150°C 
with lower values related to wind driven fires (Ventura et al., 1998). 
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The report provided by Ventura et al., (1998) presents data derived from an experiment that was carried out, 
and the report of Mendes-Lopes (1998) provides data about the same experiment. These measurements were 
made in a wind tunnel where wind velocity, fuel moisture and slope were used as parameters to study fire 
propagation in beds of Pinus pinaster needles. The experiments were carried out for 2 different fuel moisture 
contents (10 % and 18 %), 4 different slopes (0°, 5°, 10° and 15°) and 7 different wind conditions (from -3 
m/s to +3 m/s in steps of 1 m/s). It was very interesting to read the comments provided in the report on the 
relationship between temperature distribution and flame characteristics. 


The results obtained in experimental work by Ventura (1998) are plotted in Figure 1.14., where the maximum 
temperature reached and the corresponding mean spread rate for all the wind-driven situations are illustrated. It 
can be seen that the maximum temperature reached for the drier fuel is not dependent on the spread rate (but 
the rate of spread depends on the temperature). However, this is not the case for the 18% fuel when the spread 
rate is low, with an initial increase in maximum temperature for an increase in rate of spread, and levelling off 
for values similar to those of the 10% fuel at high rates of spread. Ventura (1998) considers that as the flame 
leans further forward, the amount of radiative energy that reaches the virgin fuel increases. 
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Figure 1.14 Maximum temperature and rate of spread for wind-driven flames: a) fuel moisture content: 
10%; b) fuel moisture content: 18%. (Ventura, 1998) 


Marcelli (2002 and 2004) reports experiments, which have been carried out in order to compensate the 
temperature of the flame region for a fire spreading across a pine needle bed and, hence, to measure the 


upward gas velocity (Table 1.10) in this same region. Experimental fires were conducted in a closed room 


without any air motion. The properties of the ambient air were measured for each test; the temperature 
ranged between 25 to 30°C and the relative humidity varied from 40 to 50%. The fuel was dead needles of 
Pinus pinaster. Their ratio of area to volume and density is 4550 m" and 680 kg/m’, respectively. The 
depth of fuel was approximately 2.3 cm, giving a bulk density of 21.74 kg/m’ and a packing ratio of 0.032. 
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Probe interval [m] Transit time [s] Velocity [m/s] Relative error [%] on 
estimated velocity 
0.04 — 0.06 0.065 0.306 2.04% 
0.06 — 0.08 0.027 0.747 7.85% 
0.08 — 0.12 0.043 0.926 3.28% 
0.12 -—0.16 0.034 1.160 4.40% 
0.16 — 0.20 0.030 1.321 6.80% 
0.20 — 0.24 0.030 1.342 6.18% 
0.24 — 0.28 0.026 1.532 9.34% 


Table 1.10 Mean vertical gas velocities (Marcelli et al., 2002) 


The experiments used a fuel with a fixed fuel moisture content, thus the needles had to be first dried in an 
oven at 60°C for 24 hours. They were not completely dehydrated at every experiment since they partially 
rehydrated during their removal from the oven and the ignition. Thus a sample of needles was taken for 
every run and loaded to evaluate the moisture content of the fuel before the burning. It was found that the 
needles were conditioned to a moisture content of 1% - 3%. 


Some of the results of the experimental work done by Marcelli (2002 and 2004) are plotted in Fig. 1.15 
and Fig. 1.16 


1000 


1.6 . 
© | ) 
900} ie} = 14h 
| Z fe) o) 
— | 
qe 800F - 1°) £1.25 
OS 700; 2 | 
g ° 2 ° 
= 600} = 0.8) 
e | O S fe) 
500+ 0.6 
O° 
400+ 0.4 
fe) 
‘ ‘ O., 
300 7 _ 02 
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.05 0.1 0.15 0.2 0.25 03 
Height (m) Height (m) 
Fig.1.15 Mean temperature profile Figure 1.16 Mean gas velocity profile 
(Marcelli et al., 2002) (Marcelli et al., 2002 


Download free ebooks at bookboon.com 


57 


Please click the advert 


Biomass Properties and Prediction Tools for 
Biomass Properties and Fire Prediction Tools Vegetation Wildland-Urban Interface (WUD) Fires 


Available data for temperature and velocity distribution for and above a pine needle vegetation bed 
(Ventura, 1998; Marcelli, 2002; Marceli, 2004) indicate that the maximum value for temperature during a 
fire could reach 1000°C (Fig. 1.15 and 1.16). This value is measured immediately above the vegetation 
bed. As the distance (height) increases, the temperature decreases and this corresponds to a reduction in 
flame intensity. Variations in vertical velocities feature the opposite effect with velocity increasing when 
height is increased (Fig. 1.11). From the results presented (Table 1.10), it is obvious that the vertical 
velocity of combustion gases in the vegetation bed is around 0.3 m/sec and rises up to 1.5 m/sec at heights 
of 0.25 meters. 
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